June, 1996 



On 



OO 



X 

5-H 



Effects of QCD Resummation on Distributions of 
Top— Antitop Quark Pairs Produced at the Tevatron 



g i S. Mrenna W-Q and C.-P. Yuan Wfl 



a) High Energy Physics Division, Argonne National Laboratory 
Argonne, IL 60439, U.S.A. 
(b) Department of Physics and Astronomy, Michigan State University 



Argonne, IL 60439, U.S. A 

m 

CO 

■ East Lansing, MI 48824, U.S. A 

Q^h! Abstract 

We study the kinematic distributions of top-antitop quark (tt) pairs produced 
^SZj ■ at the Tevatron, including the effects of initial state and final state multiple 

soft gluon emission, using the Collins-Soper-Sterman resummation formalism. 
The resummed results are compared with those predicted by the showering 
event generator PYTHIA for various distributions involving the tt pair and the 
individual t or t. The comparison between the experimental and predicted 
distributions will be a strong test of our understanding and application of per- 
turbative QCD. Our results indicate that the showering event generators do not 
produce enough radiation. We reweight the PYTHIA distributions to agree with 
our resummed calculation, then use the reweighted events to better estimate 
the true hadronic activity in ti production at hadron colliders. 



1 mrenna@hep . anl . gov 
2 yuan@msupa. pa. msu.edu 



1 Introduction 



Because the top quark mass m t is comparable in magnitude to the vacuum expec- 
tation value t>=246 GeV [|TJ, studying the interactions of the top quark may provide 
information on the mechanism of electroweak symmetry breaking or the generation 
of fermion masses |§. To observe any new physics effect in the top quark system, 
one has to know first the Standard Model prediction for the production rate and the 
kinematics of the top quarks produced at colliders. We concentrate on top quark 
pairs produced in hadron collisions. The next-to-leading-order (NLO) prediction for 
the production rate of ti pairs has been known for several years 0]. Since then, 
several studies [7j have extended this result to include the effect of soft gluon 

radiation on the production rate of ti pairs at hadron colliders. The NLO prediction 
of the ti production rate varies from the leading order (LO) prediction by 15-35% at 
the Tevatron for m t = 175 GeV, depending on the choice of scale for the hard scat- 
tering process, and is less sensitive (~ 10%) to the scale choice. The resummation 
of multiple soft gluon emission may increase the production rate by another 10%, 
depending on the prescription for performing the resummation || [Tj. Besides testing 
the production rate, it is also important to study the kinematics of the top quark to 
probe possible new physics associated with its production or decay. 

It is well established that the transverse momentum distribution of the ti pair 
cannot be described by the NLO perturbative calculation for small Qt- The same 
is true for the NLO prediction of the transverse momentum of electroweak gauge 
bosons M. This implies that the transverse momentum p l T of the top quark cannot 
be accurately predicted by the NLO calculation, especially for ti pairs with small Qt, 
where the data dominate. The effects of the initial state and the final state multiple 
soft gluon emission must be resummed to predict the kinematic distributions of the 
top quarks produced in ti events at hadron colliders. This work expands upon an 
earlier study of the kinematics of heavy quark pairs 11 using the Collins-Soper- 
Sterman formalism to perform the resummation |]I0" |. We closely follow the notation 
used in Ref. 0. 
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Our present understanding of the tt pair kinematics is based on showering event 
generators, such as HERWIG, ISA JET and PYTHIA [12]. A further goal is to quantify 



the successes and limitations of such generators and to make progress towards a more 
complete description of the tt pair and individual t and i kinematics This will 
have important implications for the precision measurement of the top quark mass. 

Following this introduction, we have organized this study into five additional sec- 
tions. Sec. 2 contains a review of the Collins-Soper-Sterman (CSS) resummation 
formalism. In Sec. 3, we present our numerical results for the qq — > tt and gg — > tt 
subprocesses using this formalism. We compare our results with the showering event 
generator PYTHIA in Sec. 4. Based on the results of Sec. 4, an improved estimate of 
the hadronic activity in tt events is presented in Sec. 5. Finally, Sec. 6 contains our 
conclusions. 

2 The CSS Resummation Formalism 

Soft gluon resummation has been applied successfully to predict the rate and 
kinematics of electroweak gauge boson production at hadron colliders f| HU- 
Although the applicability of the CSS formalism to tt pair production (which is a 
colored final state) has not been proven in the literature, the large top quark mass 
relative to Aqcd should suppress contributions from color configurations not included 
in the CSS resummation formalism. This should be more correct when the t (and 
t) are produced in the central rapidity region. All the leading and sub-leading loga- 
rithmic singularities associated with the initial state radiation in the NLO expression 
for qq —>■ tt production are universal to those found for electroweak gauge boson 
production M. For gg — > tt production, they are the same as those for Higgs boson 
production [J15, Obviously, there are also singularities associated with the final 



state radiation in tt production which are absent in either electroweak gauge boson 
or Higgs production. 

Our starting point for applying the CSS formalism to tt production at hadron 



3 



colliders is the resummed expression for the differential cross section: 



da{h x h 2 -> tt + X) 
^ dQ 2 dy &Q\ d4> t t d cos 9 1 



7T 



36SQ 5 



{(2^ / rf26e ^ T ''E^( & *'Q> a; i^2,^0,C' 1 ,C' 2 ,C'3)F^ p (6,Q,x 1 ,x 2 ] 



+ Y(Q T ,Q,x 1 ,x 2 ,9, ( p,C 4 )y 



In this expression, the production rate is described in terms of the mass Q, rapidity 
y, transverse momentum Qt, and azimuthal angle <j) t t of the tt pair in the laboratory 
frame, and the polar angle 9 and azimuthal angle in a special center-of-mass 
frame for the tt pair, the Collins-Soper frame |I7J . The center-of-mass energy 
hadrons h\ and h 2 fixes the parton momentum fractions X\ = x 2 = -^e~ y . The 

renormalization group invariant Wjk is given by 

W jk (b, Q, xi, x 2 , 9, 0, d, C 2 , C 3 ) = exp {-5(6, Q, C 1; C 2 )} 

(pji <g> (Xi) (C fem <g> / TO /ft 2 ) (x 2 ) + (Ch ® (Xi) ((7 im (8) /m/fta) (x 2 ) 



X 



X 



«s(C 2 Q) 
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H jk {Q, cos 9, m t ) 



(2) 



where a s is the strong coupling constant, /3i = ^(33 — 2n/) (n/ is the number of light 
quark flavors) and Cg> denotes the convolution integral 

(Cfl ® /,/ hl ) (xx) = £ |i 6, a* = y , Cx, C a )/ I/fcl (Ci, A* = y )■ (3) 

Because there are two separate hard processes in the LO calculation, there are two 
functions W qq - and W gg . The dummy indices / and m are meant to sum over quarks 
and anti-quarks or gluons, and summation on double indices is implied. The angular 
function TCjk(Q, cos 9, m t ) in Eq. @ for j = q, k = q is 



2 - (3 2 + (3 2 cos 2 9 



and for j = k = g is 

3(7 + 9(3 2 cos 2 i 



32 (1 - /3 2 cos 2 fl) 2 



1 + 2(3 2 - 2/3 4 - 2(3 2 {l - (3 2 ) cos 2 9-^ cos 4 9 



4 



where (3 = Jl — 4mf/Q 2 . We generically refer to Wjk as the CSS piece. The Sudakov 
form factor S(b, Q, C\, C 2 ) is defined as 



S(b,Q,C u C 2 ) 



ctQ 2 dp? 



r,2 



\n\^-)A(a a (fi)) + B(a a {fi)) 



P 2 



The functions A, B and the Wilson coefficients Cji were given in Ref. [11]] for qq 



(4) 



tt 



(see Eqs. (3.19) to (3.26) for B™, A^ 2 \ Cf k \ C$, and CfJ\ and in Ref. |TJ for 



^ -> tt (see Eqs. (3.8) to (3.9) for B^, Cg>, Cj', and C«). 3 In those results, 
the constants Ci, C*2 and C3 = /x6 were introduced when solving the renormalization 
group equation for the CSS piece Wjk- The canonical choice of these renormalization 
C3 = 2e~ 7E = 60 and C 2 — 1 [[11], [l5| . (7s is the Euler con- 



constants is C\ 

stant.) To test the dependence of our numerical results on the particular choice of 
the renormalization constants, we consider the set of constants such that C\ = C 2 b 
and C3 = bo. This choice eliminates large constant factors in the expressions for the 
A, B and Cjk functions. Because the tt final state is colored, there is an additional 
contribution to the B function inside the Sudakov factor due to final state gluon ra- 
diation. The mass of the top quark regulates a potential collinear singularity so that 
the final state contributes only ln^(^r) terms due to soft gluon emission. Since there 
are no ln^ 

contribution to the B function was given originally in Ref 

(i) 



r) contributions at NLO, there is no A^ nal function. The additional 



B 



final 



c F 



(5) 



V+P, 

Near threshold, when p -> 0, B)^ nal = -C F with C F = 4/3 in QCD. 

As shown in Eq. @, the upper limit of the integral for calculating the Sudakov 
factor is p, = C 2 Q, which sets the scale of the hard scattering process when evaluating 
the renormalization group invariant quantity Wjk, as defined in Eq. (Q). The lower 
limit JX = Ci/b — bo/b determines the onset of non-perturbative physics. 
The F-term in Eq. ([I]) is defined as 



Y(Q T , Q, x lt x 2 , 9, 0, C 4 ) = f 1 ^ I" ^ f) 

Jxi qi Jx 2 q, 2 N=1 



u s {C A Q) 



7T 



(2+N) 



3 The superscripts (0), (1), and (2) represent the order in a s . A,B and Cjk are all calculated in 
the MS (modified minimal subtraction) scheme. 



x/ Vhl (£i; C,Q) rZ\Qt, Q, ^, ^, 5, 0) / m/ta (6; C 4 Q), (6) 

where the functions only contain contributions less singular than ^4- x (1 or ln(^-)) 
as Qt — > 0. We denote those singular contributions as the singular-piece in contrast 
to the regular V-piece. The scale of the F-piece is specified by the choice of C4. To 
optimize the perturbative expansion, that is, to minimize the contribution of logarith- 
mic terms ln(C4) from higher order corrections, we choose C4 = 1 in calculating the 
F-piece. More specifically, to obtain the regular F-piece, we subtract the singular- 
piece for qq — > tig, gq/q — > tiq/q, and gg — > tig (which can be obtained by expanding 
Eq. 1 to order with C\ = C2&0 and retaining those terms proportional to Q? 2 ) 
from the squared amplitude for the tree level processes qq — > tig, gq/q — > tiq/q, and 
99 -> tig. 

In Eq. ([!]), the impact parameter b is to be integrated from to 00. However, for 
b > b max , which corresponds to an energy scale less than l/b max , the QCD coupling 
a s becomes so large that a perturbative calculation is no longer reliable. 4 The non- 
perturbative function F NP is needed in the formalism with the general structure 



Fjk (b,Q,Q ,xx,x 2 ) = exp 



(Q 2 \ 

-In — 2 ]h 1 {b) - h j/hl {xi,b) - h k/h2 (x 2 ,b) 



(7) 



The functions h±, hj/^ and h k / h2 cannot be calculated using perturbation theory and 
must be measured experimentally. Furthermore, the CSS piece W is evaluated at 
with 

b 



;i + (P/b max y 

such that fe* never exceeds b max |J. 

To obtain the final product of our calculation, the kinematics of the t and i, we 
transform the four-momentum of t (= p M ) and i (= p M ) from the Collins-Soper frame 
to the laboratory frame. The resulting expressions are:Q 

if = 9- + sin 6 cos (frX" + sin 6 sin W + cos QZ% 



4 We use b max = 0.5 GeV 1 in our calculation. 
5 Our convention is = (q°, q 1 , q 2 , q 3 ). 



p 1 * = g M — jf, 



g M = (Mt cosh y, Qt cos (f), Qt sin 0, My sinh , 
O M 2 



' {() ■ a 1 ' () //' 



with Q± = ± g 3 ), Q = V¥, M T = ^Q 2 + Q T , y = \ ln(g±), n" = ^(1, 0, 0, 1), 

and n v = ^(1,0, 0,-1). 

3 The Numerical Results of Resummation 

In this section, we present numerical results for the qq — > tt and gg — ► sub- 
processes after applying the resummation formalism outlined in the previous section. 
For these results, we have assumed m t = 175 GeV for tt production at the Tevatron 
l a PP 

collider) with y/S = 1.8 TeV. 
As explained in the previous section, the CSS piece depends on the renormalization 
constants C\,C 2 = Ci/bo and C 3 = &o- The choice of C 2 indicates that the hard scale 
of the process is Q = C 2 M t t, where M t i is the invariant mass of the tt pair. We use 
CTEQ3M NLO parton distribution functions (PDF's) the NLO expression for 
a s , and the non-perturbative function JOS] 



F NP (b, Q, Q , x u x 2 ) = exp j-^ 2 - g 2 b 2 In ^ J - 9l g 3 b\n (100a*r 2 ) J , (9) 

where g x = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , g 3 = -1.5 GeV" 1 and Q = 1.6 GeV.Q Finally, 
the CSS piece is fixed by specifying the order in a s of the A, B and Cjk functions. 
We adopt the notation (M, N) to represent the order in a s of #),BW andCjf. 
The choice (1, 0), for example, means that A and B are calculated to order a s , while 

6 These values were fit for CTEQ2M PDF and C 2 =l, and in principle should be refit for CTEQ3M 
PDF and different values of C2. Also, for the gg process, 32 should be replaced by 32— f^- — 32 1 via 

A qq 

the renormalization group argument for the In(g^) dependence of the non-perturbative function. 
Since the gg channel is numerically less important at the Tevatron, we still use Eq. 9 in this study. 
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Cjk(z) is either or 5(1 — z) depending on j and k. Also, (1,0) means the a^Hjk 
term in Eq. is not included. 

If the CSS piece is expanded to order afi, then it contains all the ln(Sj-) terms 
predicted by the NLO calculation. As discussed in the previous section, the regular 
V-piece is calculated by subtracting all the singular contributions, which grow as 
-4r x (1 or ln(Sr)) as Qt — > 0, out of the NLO tree level processes qq — ► itg, (?g/ g — > 
tiq/q, and — > tig. In Fig. 1, we show the relative sizes of the CSS piece, the 
singular-piece expanded to order a^, and the NLO tree level result (which is also 
order afy as a function of Qt- The regular V-piece, which is defined as the difference 
between the NLO tree level result and the singular-piece, is small for Q T up to 50 
GeV. Its relative contribution starts at and reaches about 30% at Qt=50 GeV. The 
two curves with singular behavior as Qt — > have been cut off at Qt="2 GeV for the 
purposes of this figure only. 

We conclude that the F-piece is not important for small Qt up to about 25 GeV, 
while most of the rate occurs at much smaller Qt- In Table p], we present the rate 
for the CSS piece alone for Qt < 50 GeV for each choice of order (M, N) and the 
dependence of this rate on the renormalization constant C2. Results for the gg — > ti 
channel are only presented up to order (1,1). In the same table, we show the mean 
and standard deviation for each order (M,N). For the highest order, (2,1), the 
variation with the hard scale set by C2 of the qq — > ti channel is only 2%. For 
the gg — > ti channel, there is only a marginal improvement in the variation from 
higher order, though the overall correction at higher order is large. This raises some 
concern regarding the higher order dependence of this result. A similar behavior is 
exhibited in the gg channel for the other resummation schemes. Fortunately, at the 
Tevatron, the dominant contribution to the ti production rate (~ 90%) comes from 
the stable qq — ► ti channel. In Fig. 2, we show the relative contributions of the qq 
and gg channels to the CSS piece in this formalism. The total rate for ti production 
is obtained by adding the F-piece to the CSS piece. These results are compiled in 
Table |2] in a similar fashion as in Table [TJ We have also included a column a[ N+2 ^ a per t 



8 



which shows the LO (for N = 0) and NLO (for N = 1) perturbative cross sections for 
comparison. The choice of order (1,1) for both the qq — > ti and gg — > ti contributions 
represents the same order as the NLO calculation. With the choice of hard scale 
Q = M t t/2 (C*2 = 1/2), the integrated total production rate for top quark pairs is 
found to be 5.64 pb, which agrees within 10% with the NLO result 5.06 pb evaluated at 
the similar hard scale Q = m t . This indicates that the CSS resummation formalism 
presented in the previous section contains the dominant contribution of the NLO 
result to the production rate of ti pairs. Our average resummed result 5.58 ± 0.09 pb 
agrees with Ref. which was obtained using principal value resummation. In the 
literature, the cross section for ti production is usually given by taking the hard scale 
to be a multiple of m t rather than M#. For the LO and NLO calculations and other 
resummation formulations which implicitly integrate out the Qt dependence, the only 
hard scale left in the problem is the mass mt. The LO and NLO perturbative results 
in the final column of Table ||] are evaluated at the scale Q = Ci x (2m t ). The scale 
Q = m t , then, corresponds approximately to choosing the renormalization constant 
C*2 = 1/2 in the CSS formalism. 

Our resummation calculation only includes the finite contributions from the vir- 



tual diagrams which are the same as those in the Drell-Yan process |TT| for the 
qq — > ti channel and those in Higgs production [OJ for the gg — ► ti channel, plus 
terms containing the running of a 2 s in the hard part of the cross section multiplying 
7~Cjk- Since there is no final state QCD radiation in the Drell-Yan or Higgs produc- 
tion processes, this is clearly an approximation. If the exact virtual corrections were 
included, the Wilson coefficient functions Cjf. would have to be modified in a manner 
consistent with the CSS formulation, i.e. assuming that the initial and final state 
gluon radiation for the ti process factorizes in a similar fashion as the initial state ra- 
diation alone in the Drell-Yan process. While such a factorization is reasonable, there 
is no formal proof. Therefore, we approximate the qq —>■ ti (gg — > ti) finite virtual 
corrections with those from the Drell-Yan (Higgs production) process. Furthermore, 
the factor for qq — > ti inside the Sudakov factor is likely to be different from that 
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Table 1: CSS Contribution to the Total Cross Section for Top-Antitop Production 
at the Tevatron 
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Table 2: Total Cross Section for Top-Antitop Production at the Tevatron in the CSS 
Formalism 



10 



for the Drell-Yan process, so we did not include it in our calculations. As shown in 
Tables [l] and 0, as C2 varies between 1/4 and 1, the total rate varies by only a few 
percent, which is about the same magnitude as the uncertainty in the NLO calcula- 
tions. Clearly, all the other finite contributions from the virtual diagrams other than 



those similar to the Drell-Yan virtual contributions for qq — > tt are small ||22j| . Since 
for this process is the same as for the Drell-Yan process, we have included it in 
our order (2,1) calculation for qq — > tt to improve our predictions, which yields the 
results shown in the first row of Table 2. 

Because of the agreement of our predictions for the total rate with the NLO cal- 
culation, we apply the results of our calculations to study the kinematic distributions 
of the tt pair and the individual t or t produced in hadron collisions. In the next sec- 
tion, we present these kinematic distributions and compare them to those predicted 
by showering event generators. We wish to stress that this is the only type of compar- 
ison that is sensible, because the LO calculation predicts a 5{Qt) dependence for the 
Qt distribution, the NLO calculation cannot accurately describe the Qt <C m t region 
of phase space, and other resummation formalisms integrate out the Qt dependence 
and, thus, cannot predict the kinematic distributions of t and t. 

In the following sections, we present results only for order (2, 1) and the canonical 
choice of renormalization constants C%—\. There are at least two reasons for doing 
this. First, the coefficients of the non-perturbative function used in this study were 
fit to data assuming C% — 1. While the non-perturbative function can affect the 
shape of distributions, it does not affect the total rate. Therefore, the stability of 
our results to variation of C2 is still valid, but we cannot trust other choices of C2 
to give the correct shape. Second, since we do not integrate out the kinematics, we 
cannot argue that the only scale left in the problem is m t . Since an s-channel process 
dominates, Q = M t i is motivated by the dynamics. 

We have also checked the effect of neglecting final state radiation, and have found 
that it reduces the total contribution of the CSS piece to tt production by about 10% 
in our approximation and slightly changes the shape near the peak. Without the final 
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state radiation, however, the F-piece is not finite as Qt — > 0; therefore we include it 
in our results. 

4 Comparison with the Showering Monte Carlo 
Technique 

In this section, we compare our resummed results for the kinematics of the tt 
pair and the individual t or i with those predicted by the showering event generator 
PYTHIA. To make this comparison more understandable, we first explain the approx- 
imate implementation of resummation in such generators. The starting point for the 
showering Monte Carlo technique is the observation that the leading logarithmic sin- 
gularities in NLO calculations are contained in the Altarelli-Parisi splitting kernels. 
In this leading log approximation, successive parton emissions occur independently, 
modulo some angular ordering effects. In a Monte Carlo simulation, which has ex- 
plicit finite cutoffs for the energy of the emitted radiation, it is possible to treat these 
emissions as a Markov chain stretching from the hard scattering backwards to the 
initial state partons ||20|| . Essentially, one chooses the kinematics for a process at the 
hard scale Q max based on the LO cross section, then calculates the probability for no 
parton emission in evolving from a high scale t max rs ln(<5^ aa ./ Aq CD ) to a lower scale 
t. This probability is given in the leading log approximation by the Sudakov form 
factor e ~ s ( x ' t,tmax \ where 

rt max r q s ({) xf a (x,t) 

= dt}_^ dz— — -^-P a ^ bc {z). 

Jt a J 27T xf b (X,t) 

Here, f a {x,t) and fb{x,t) are parton distribution functions for partons a and b and 
Pa~,bc{z) is the Altarelli-Parisi splitting function for the branching a — > be with 
momentum fraction z. In improved treatments of the Sudakov form factor, such as 
in PYTHIA, a s is evaluated not at the scale Q 2 , but at Q 2 (l — z) J2TJ. If no radiation 
occurs down to some cutoff t m i n , ~ 2 GeV in PYTHIA, then the parton is placed on the 
mass shell. On the other hand, if it is determined that radiation does occur before the 
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cutoff, then a new parton is emitted. The type of branching a — > be which led to the 
parton emission is determined by the relative weights of the x integrals over P a ^bc- 
The initiator of the splitting is given the virtuality t, and the process continues until a 
parton reaches the scale t m j„, and it is placed on the mass shell. With each splitting, 
the remaining kinematics are sampled so that energy and momentum are conserved. 
The result is a total cross section given by the LO calculation, but with the kinematic 
distributions of a resummed calculation. In this sense, it corresponds to a choice of 
(M, N) = (1, 0) in the CSS formalism.[] Final state radiation is implemented in a 
similar fashion, but it occurs forward from the hard scattering process to the final 
state and is not weighted by the parton distribution functions. 

First, we present a comparison of the resummed and showering Monte Carlo kine- 
matics for the ti pair. Several observables which in principle cannot be extracted 
from a NLO calculation are the distributions of the transverse momentum of the ti 
pair Qt, the opening angle between t and i in the azimuthal plane A0 t j, and the 
variable z = — ^f^pi'^)^ ■ The resummed results are separately shown as solid 
lines. Displayed on the same plots are the shapes predicted by PYTHIA for two choices 
of hard scale, Q 2 = s and Q 2 = \Jm 2 + p T , where px is the transverse momentum 
of the top quark at LO. As illustrated in Table 2, the LO rate is highly sensitive to 
the choice of hard scale. The resummed estimate of the total rate is more reliable. 
Therefore, to compare the shapes of distributions, we have renormalized the PYTHIA 
results to have the same total rate as the resummed calculation. 

Fig. 3 shows the distribution of the resummed variable Qt- Note that the re- 
summed distribution is significantly harder than the showering Monte Carlo result, 
implying that the resummed calculation predicts more overall hard radiation. This 
is further demonstrated by the A0 t j distribution, the difference in azimuthal angle 
between t and i, in Fig. 4, which is depleted near A0 = tt in comparison to the 
showering Monte Carlo. Because of initial or final state radiation, the t and i are 
not expected to be exactly back-to-back. The z distribution, shown in Fig. 5, is 

7 The reason there is not exact agreement between the LO and (1,0) rates in Table 2 is because 
of the explicit cutoff b max in the Fourier transform. 
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also shifted away from a back-to-back configuration (z — 1) in the resummed result. 
Finally, the distributions of the rapidity of the it pair y t i and the invariant mass M t i 
are shown in Figs. 6 and 7 respectively. The rapidity distribution is more central 
for the resummed result and the invariant mass favors higher masses and is broader. 
These differences arise because PYTHIA only contains LO matrix elements, while the 
resummed result contains the dominant piece of the NLO correction. 

Second, we present a comparison of the kinematics of the individual t or t. Fig. 8 
shows the distribution of the transverse momentum of the individual t p T , while Fig. 
9 shows the difference in their rapidity Ay*, where Ay* = y* — y*. From the LO 
calculation, we know the scale of p T is set by m t (p T ~ m t /3 ~ 60 GeV), while 
the typical transverse momentum Q T is much smaller. Therefore, there is not much 
difference in these distributions. Likewise, Ay* is more sensitive to the PDF (which 
determines the boost of the ti pair) than the transverse momentum Qt, so we do 
not expect to observe a large difference. In conclusion, the showering generators (as 
typified by PYTHIA in our study) reproduce the CSS distributions for the individual 
t and i and the tt pair kinematics in our plots to a 10% level bin-by-bin, although 
the overall shapes are generally different and the complete resummed results indicate 
more overall hard radiation. Based on these results, we seek to improve the showering 
Monte Carlo technique by using our knowledge of the resummed Qt distribution. This 
is discussed in the next section. 

5 Jet Activity in ti Events at Hadron Colliders 

Despite its limitations in predicting the correct rate, the showing Monte Carlo 
technique has a mechanism to approximate the complete resummed result for the Qt 
distribution. Furthermore, the showering event generator gives a phenomenologically 
accurate description of all the details of the event. The resummation calculation only 
predicts the vector sum of all soft gluon radiation, but has no power to predict how 
the radiation is distributed amongst individual gluons or quarks. Such details are 
crucial for estimating the amount of jet activity, which will affect the determination 
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of the top quark mass by reconstructing jets from top quark decays as done by CDF 
and DO at the Tevatron These estimates are used to engineer cuts to enhance 
the signal and to optimize the choice of jet definition. In this section, we present a 
simple synthesis of the showering Monte Carlo technique and the full resummation 
calculation to realize a better estimate of the hadronic activity This synthesis is 
accomplished by reweighting events generated by PYTHIA to agree with the resummed 
Qt distribution and rate. 

To demonstrate the improved predictive ability of the resummed approach, we 
present the distributions of the jet multiplicity and scalar sum of the transverse 
energy in ti events. We have isolated the contributions from initial and final state 
radiation only, so that none of the the t or i decay products (or additional QCD 
radiation contributions from them) are included in these plots. We define jets using 
a simple calorimeter simulation (with cell segmentation At] x A0 = .1 x .1) and a 
clustering algorithm based on a cone size R = .4 and a minimum jet transverse energy 
EJp m =5 GeV. The energy E deposited in each cell of the calorimeter is smeared 
with a Gaussian resolution ^ = -^=. The jet multiplicity distribution is shown in 
Fig. 10. Note the shift in the peak value of the distribution from jets to 1 jet. The 
resummation based Monte Carlo clearly predicts more hard radiation in ti events 
at hadron colliders. Similar information is conveyed in Fig. 11, which shows the 
scalar sum of the transverse energy of the jets defined as above. While the details of 
jet observables can only be studied using the showering Monte Carlo technique, our 
knowledge of the resummed Qt distribution allows improvement. It is straightforward 
to extend these results to include the decay products of the t and i and even the hard 
gluons from the QCD radiative decay of the top quark |23|. We shall leave this for 
further study. 

As a final point of comparison, we have attempted to quantify the effect of hard 
gluon radiation on the extraction of m t from data. The issue at hand is how often 
a hard gluon from radiation is misidentified as a top quark decay product. Using 
a sample of events where t(— > bW + (— > e + v e ))i(— > bW~(-^ jj)), we cluster particles 
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into jets as described above and analyze events with 4 or more hard jets [E 3 T >15 
GeV, | r}W | < 2). We separate the jets than can be identified with a t decay product 
(for simplicity, we assume the b is correctly tagged) and find the one with the lowest 
Et- The fraction of events where another jet (i.e. one not from t decay) has a higher 
E T than this is an estimate of the importance of the hard radiation. While this is only 
a crude estimate of the real effect, we find that the fraction (= .56) does not differ 
between the standard and improved PYTHIA. This implies that the hard gluon error 
on the top mass measurement is well estimated by showering generators, though this 
requires further study. We have not attempted to quantify the effect of soft gluon 
radiation, where the radiation is not resolved as an individual jet but can overlap 
with the ti decay products, though the complete resummed result indicates that this 
too should be enhanced. 

6 Conclusions 

To further study the interactions of the top quark and to better measure its 
mass, we must understand the kinematic distributions of the transverse momentum, 
rapidity and azimuthal angle of the top quarks. The kinematics of the top quarks 
produced at hadron colliders can be accurately predicted only after resumming the 
multiple soft gluon emissions in either the initial state or the final state. In this work, 
we have adopted the CSS resummation formalism to obtain the kinematics of the 
top quarks. The approximation we made in this study should be adequate because 
of the large top quark mass. The important consequence of the large top quark 
mass is that the logarithmic terms from the initial state are more important than 
those from the final state. In the former case, there is double-logarithmic behavior 
due to both soft and collinear gluon radiation, while, in the latter case, only single- 
logarithmic behavior due to soft but not collinear gluon radiation is possible up to 
the order a^. We compared our resummation calculation with those predicted by 
the full event generator PYTHIA and found that the latter does not give the same 
prediction as ours. In view of the fact that the full event generator has been widely 



16 



used by our experimentalist colleagues for analyzing their data, it is important to 
point out the difference between the results from this approach and those from the 
analytical calculation. In determining the mass of the top quark, for example, one 
needs to account for the jet activity in at least two cases: (1) when a hard jet from 
the showering is misidentified as one of the top decay jets, and (2) when the soft 
radiation is included in the energy determination of true top decay products. The 
hybrid approach presented in Section 5, which relates the transverse momentum of the 
top-antitop quark pair from an analytic calculation to the balancing gluon radiation 
from a showering Monte Carlo, contains a more realistic description of event structure, 
which could be used for choosing kinematic cuts and tuning the jet energy correction 
algorithm. 

Acknowledgments 

We thank Ed Berger for starting us on this project and providing helpful com- 
ments. C.-P. Y. also thanks the CTEQ collaboration and C. Balazs for discussions. 
This work was supported in part by NSF grant PHY-9309902 and by DOE grant 
DE-FG03-92-ER4070 1 . 

References 

[1] F. Abe et al, Phys. Rev. Lett. 73, 225 (1994); 
S. Abachi et al, Phys. Rev. Lett. 72, 2138 (1994). 

[2] Ehab Malkawi and C.-P. Yuan, Phys. Rev. D50, 4462 (1994). 

[3] R.S. Chivukula, E. Gates, E.H. Simmons and J. Terning, Phys. Lett. B311, 157 
(1993); 

[4] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303, 607 (1988); B327, 49 
(1989); 

W. Beenakker, H. Kuijf, W.L. van Neerven and J. Smith, Phys. Rev. D40, 54 



17 



(1989) ; 

R. Meng, G.A. Schuler, J. Smith and W.L. van Neerven, Nucl. Phys. B339, 325 

(1990) . 

[5] E. Laenen, J. Smith and W.L. van Neerven, Nucl. Phys. B369, 543 (1992). 

[6] E.L. Berger and H. Contopanagos, Phys. Lett. B361, 115 (1995). 

[7] S. Catani, M.L. Mangano, P. Nason and L. Trentadue, CERN-TH/96-21, Jan- 
uary 1996. 

[8] G. Altarelli, R.K. Ellis, M. Greco, G. Martinelli, Nucl. Phys. B246, 12 (1984). 

[9] E.L. Berger and R. Meng, Phys. Rev. D49, 3248 (1994). 

[10] J. Collins and D. Soper, Nucl. Phys. B193 381, (1981); Erratum B213 (1983) 
545; B197, 446 (1982). 

[11] J. Collins, D. Soper and G. Sterman, Nucl. Phys. B250, 199 (1985). 

[12] G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour, and 
L. Stanco, Comp. Phys. Commun. 67, 465 (1992); H. Baer, F. Paige, S. Pro- 
topopescu, and X. Tata, "Simulating Supersymmetry with ISAJET 7.00," FSU- 
HEP-930329, SSCL-Preprint-441 (1993); T. Sj6strand, Comp. Phys. Com- 
mun. 82, 74 (1994). 

[13] S. Frixione, M.L. Mangano, P. Nason, and G. Ridolfi, Phys. Lett. B351, 555 
(1995). 

[14] PB. Arnold and R.P Kauffman, Nucl. Phys. B349, 381 (1991). 

[15] C.-P. Yuan, Phys. Lett. B283, 395 (1992). 

[16] R.P. Kauffman, Phys. Rev. D45, 1512 (1992). 

[17] J. Collins and D. Soper, Phys. Rev. D16, 2219 (1977). 



18 



[18] H.L. Lai, J. Botts, J. Huston, J.G. Morfin, J.F. Owens, J.W. Qiu, W.K. Tung, 
H. Weerts, MSU preprint MSU-HEP-41024, Oct 1994. 

[19] G.A. Ladinsky and C.-P. Yuan, Phys. Rev. D50, 4239 (1994). 

[20] T. SjSstrand, Phys. Lett. 157B, 321 (1985). 

[21] A. Bassetto, M. Ciafaloni and G. Marchesini, Phys. Rep. 100, 202 (1983). 

[22] R. Meng, G.A. Schuler, J. Smith, and W.L. van Neerven, Phys. Lett. B339, 325 
(1990). 

[23] S. Mrenna and C.-P. Yuan, Phys. Rev. D46, 1007 (1992). 



19 




Figure 1: The relative importance of the F-piece with respect to the CSS piece as a 
function of Qt- 
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Figure 2: The relative contributions of qq — > tt and gg — > tt to the CSS piece as a 
function of Qt- 
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Figure 3: The transverse momentum of the tt pair for C2—I and for different 
choices of hard scale for the showering Monte Carlo PYTHIA. The PYTHIA rate has 
been renormalized to the CSS rate. 
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Figure 4: The difference in azimuthal angle between the t and t A0 t f. 
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Figure 5: z distribution of the t and t. 




Figure 6: Rapidity of the tt pair y t i- 
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Figure 9: Rapidity difference between the t and t Ay t . 
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Figure 10: Jet Multiplicity from initial and final state radiation as predicted by 
PYTHIA (dashed line) and the hybrid PYTHIA- CSS resummation (solid line). 
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Figure 11: Scalar sum of jet transverse energy Et from initial and final state radia- 
tion as predicted by PYTHIA (dashed line) and the hybrid PYTHIA- CSS resummation 
(solid line). 
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